function [ hasil ] = y_Adams_Bashforth_Moulton( f, x0, y0, b, h )
%Y_ADAMS_BASHFORTH_MOULTON Summary of this function goes here
%   Detailed explanation goes here
    n = (b - x0) / h;
    y = y0;
    
    y1 = y_RK3(f, x0, y0, x0 + h, h);
    y2 = y_RK3(f, x0, y0, x0 + 2 * h, h);
    y3 = y_RK3(f, x0, y0, x0 + 3 * h, h);
    x = x0 + 3 * h;
    
    syms a b;
    
    for r = 4 : n
        y = y3 + h / 24 * (-9 * subs(f, [a b], [(x-3*h) y0]) + 37 * subs(f, [a b], [(x-2*h) y1]) - 59 * subs(f, [a b], [(x-h) y2]) + 55 * subs(f, [a b], [x y3]));
        y = y3 + h / 24 * (subs(f, [a b], [(x-2*h) y1]) - 5 * subs(f, [a b], [(x-h) y2]) + 19 * subs(f, [a b], [x y3]) + 9 * subs(f, [a b], [(x+h) y]));
        y0 = y1;
        y1 = y2;
        y2 = y3;
        y3 = y;
        x = x + h;
    end
    
    hasil = y;
end

